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Abstract 

We prove the existence and uniqueness of solution to a classical creep dam- 
age problem. We formulate a sufficient condition for the problem to have a 
unique smooth solution, locally in time. This condition is stated in terms of 
smoothness of given data, such as solid geometry, boundary conditions, ap- 
plied loads, and initial conditions. Counterexamples with an arbitrary small 
lifetime of a structure are also given, showing the mechanical interpretation of 
imposed smoothness conditions. The proposed theory gives a rigorous frame- 
work for a strain localization analysis. The influence of the damage gradient 
on the strain localization process is characterized within this framework and 
a measure of the damage localization is proposed. 
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1 Introduction 

Structures made of metals and alloys are often used in different brunches of industry 
at elevated temperatures (higher than 0.3 times the melting temperature). Typical 
examples are pressurized pipes and vessels in power and chemical plants, gas turbines 
and so on. Even subjected to moderate loads these structures experience irreversible 
creep deformations which influence the stress response in long time scales. The 
lifetime of such structures is limited by damage processes induced by the nucleation 
and growth of microscopic cracks and cavities. The finite element method (FEM) 
is commonly used for numerical analysis of nonlinear creep-damage response in the 
framework of continuum damage mechanics (see ^0]) to estimate the remaining 
lifetime. The classical creep models are described and analyzed in the following 
monographs [HI, [1^1, dl, EH- 

The creep behavior is divided into three stages. The initial stage is characterized 
by hardening behavior with decreasing creep strain rate. The second stage is the 



stationary creep with a constant creep strain rate. The last stage is the tertiary 
creep characterized by increasing creep strain rate and a dominant softening of the 
material followed by a complete rupture. The most popular constitutive law for the 
second stage was proposed by Norton [22j and postulates the stationary creep rate as 
a power law function of the stress tensor. This constitutive law is modified by use of 
time or hardening parameters (see, for example, ^7\, ^Tj) to take the primary creep 
into account. A new internal continuity parameter i|) was introduced in the original 
work of Kachanov ^3] to simulate the material damage within the tertiary creep. 
This continuity parameter is often replaced by a dual variable, namely, Rabotnov's 
damage parameter cu = 1 — 1|> (see |13j). Within Kachanov- Rabotnov's approach 
the damage rate is postulated as a function of the stress, the temperature and 
the current damage state. This is regarded as a foundation of continuum damage 
mechanics (CDM). 

One of unsolved problems of computational CDM is the spurious mesh- 
dependence of FEM simulations (compare [TH], [20], |21]; 123) which leads to 
physically unrealistic results. Therefore numerous regularization techniques were 
proposed to prevent this mesh-dependence (see, for example, pOj). 

Proposed material models and regularization techniques are generally tested by 
series of numerical experiments. At the same time the mathematical treatment of 
nonlinear material models is very poor. Some mathematical results are given in 
|2]; [ni; |2I]- In [2] the local existence and uniqueness is proved for a coupled 
creep-damage model assuming the elastic properties are not influenced by damage 
evolution (partly coupled approach). 

For most of the used damage models it is not clear whether the corresponding 
boundary value problems are well-posed. We say that a given problem is well-posed 
(see, for example, 0) if the problem in fact has a solution; this solution is unique; 
and the solution depends continuously on the data given in the problem. In case of 
a creep-damage problem such given data are solid geometry, boundary conditions, 
applied loads, initial conditions, and material constants. A mathematically con- 
sistent problem statement is necessary for justification of analytical (see the paper 
[20]) and numerical techniques. Particulary it specifies how the difference between 
exact and approximate solution can be measured and what kind of perturbations of 
given data are allowed. 

The proper mathematical analysis of nonlinear damage models is complicated 
by instabilities due to loss of ellipticity of the corresponding differential operator 
(compare [12], [H], [2] for example). On the other hand, bifurcation does not 
happen before the appearance of completely damaged zone with cu = cu*, where cu* 
is a critical damage value (see, for example, [H])- The period of time required by the 
structure to reach this state is called crack growth initiation time t*. Therefore we 
prove existence and uniqueness of solution in sufficiently small time interval before 
failure initiation. On this time interval the deformation is stable and the problem 
can be posed correctly. Thus, the analysis of crack propagation lies beyond the 
scope of this article. 

The article is organized so that technical details of the proof do not obscure 
the main points. First, we introduce an initial boundary value problem for fully 
coupled creep-damage model. In the following section we give basic definitions of 
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function spaces, that are necessary for the formulation of the main result. In section 
4 the main existence and uniqueness theorem is formulated. One counterexample is 
provided, which illustrates the effect of damage localization. Finally we prove the 
main result and summarize our main conclusions. 

2 Constitutive equations 

Let O C be a bounded domain which represents the solid. In this work we confine 
ourselves to the plane-stress two-dimensional case. But the theory proposed here can 
be easily generalized to three dimensions. Let us assume a stationary temperature 
field. Therefore the constitutive equations do not depend on temperature. 

2.1 Fully coupled damage model 

Suppose that the damage evolution is controlled by the von Mises equivalent stress. 
Accordingly to the classical Kachanov-Rabotnov concept the constitutive equations 



for secondary and tertiary creep are summarized as follows 

(r = C"'(£-£n inax[0,T], (1) 

e^^^^As (ffvM]""\l -cu)-" inax(0,T), (2) 

cu = B (CTvm)™(1 -cu)-^ in CI X {0,1), (3) 

1 /3 

s = cr- -tr(cr)I, o-^m = \/^s:s, (4) 



where cr is the stress tensor, is the fourth-rank tensor depending on Rabotnov's 
damage parameter cu, ( ) is the time derivative, e'^'" is the creep strain, s is the stress 
deviator, o"vm is the von Mises equivalent stress, I is the second rank unit tensor, 
and A, B, n, m, q are material constants. The influence of the damage on elastic 
properties is given by the equation (see [1^], [TH] ) 

^" = 0(1-0)). (5) 

Here C denotes the tensor of linear elasticity of undamaged solid. C is linear, 
symmetric, positive definite mapping. 

Furthermore, we consider equilibrium equations 

V-cr = -q inax[0,T], (6) 

and strain-displacement relations 

£ = l(Vu + VuT), (7) 



3 



where e is the hnearized strain tensor and u is the displacement vector. The quanti- 
ties (u, e'^^, cu, cr) depend on the space variable x G H and the time variable t G [0, T] 
for some T > 0. The system is completed by boundary and initial conditions 



u = u* on 9a X [0,7], 
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lt=0 — £-0 



cr 



in CI, 



(9) 



a)|t=o = cuo 



in a. 



(10) 



Remark. The constitutive equation for damage evolution Q was generalized by 
Hayhurst (see [TT]) so that do depends on the combination ocai + (1 — (x)CTvm of 
the maximal principal stress and von Mises equivalent stress. But the proof of the 
main result {Theorem 4-i) is essentially based on the smoothness of constitutive 
equations, therefore we do not analyze this popular model here. 

2.2 Remark on the partly coupled damage model 

Equation (0) of the fully coupled model is based on Kachanov's concept of reduction 
of the effective load carrying area. This equation is a fundamental form of elastic- 
damage coupling. 

Within the partly coupled approach the influence of damage on the elastic prop- 
erties is neglected and equation (0) is replaced by 



where cu* stands for critical damage state. 

In most of engineering applications this approach is used as a simplified variant of 
the fully coupled relations in order to decrease the computational effort. Unlike the 
fully coupled model, the partly coupled approach does not require a modification 
and decomposition of the stiffness matrix on each time or iteration step. As it 
was observed in jTHI and the partly coupled approach gives a good estimation 
of failure time for some specimens and initial conditions. Nevertheless, as it will 
be shown later, this simplification should be used carefully. This model does not 
take into account the stress concentrations caused by damage inhomogeneity. For 
instance, if A = in © and B ^ in 0, then the partly coupled system (fTTj) 
describes hnear elasticity of homogeneous solid. 

Existence and uniqueness for partly coupled model were proved in ^ in case of 
thin-walled structures. Thus, the plane stress was covered as a special case of shell 
geometry. We generalize the existence proof, given in j2j, to take the fully coupled 
damage model into account. 

3 Basic notations 

The creep-damage problem (0) — (fTUI) can be formulated in a well-posed manner 
with the help of suitable function spaces. Field variables which describe the structure 




(11) 
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are considered to be elements of these infinite dimensional spaces. The corresponding 
function norms should take into account the physical essence of the problem and 
the properties of the system of equations. 

3.1 Definition of function spaces 

Let B be a Banach space endowed with a norm || ■ ||b and T be a positive real number. 
We introduce a space of continuous B-valued functions defined on the interval [0, T]. 
Definition 1 

C°([0,T],B) :={(p : [0,1] ^ B, cp is continuous}. (12) 

This space is a Banach space equipped with the norm 

||u||b,oo = sup{||u(t)||B : t G [0,T]}. (13) 

Let k G No and p ^ 1. We define the usual Sobolev space W^'^(n) (see, for 
example, [H, P, [H). 
Definition 2 

W^'P(a) := {u G Lp : G Lp for all oc = (cci, ccz) G N^, cxi + CX2 < k} (14) 
endowed with the norm 

l|ulk.p=( Y. l|D"u||9^. (15) 

aeN^, |a|<k 

Here C^i are generalized derivatives of the order |a| = ai + (Xz- 

Beside the Sobolev space W we will need a proper subspace Wo'^[D.) C W^'^p 

which is defined as follows. 
Definition 3 

Let C^(0) := {cp G C°°(n),supp cp C O} be the set of smooth functions that 
vanish near the boundary dCl. Then 

W^'^[CL):=C^{C>] (16) 

is the closure of C^(0) with respect to the norm ||u||ic,p. Note that functions from 
Wq'^(O) vanish on the boundary in the trace sense (see Definition 4)- 
Theorem 3.1 (Imbedding theorem, see Theorem 7.26 in [7]) 
Let p > 2 and D. be a Lipschitz domain in M^. Then W'^'^(n) is continuously 
imbedded m C°'^-2/P(a). 

Corollary 3.1 (Sobolev inequality) 

Let p > 2 and D. be a Lipschitz domain in M^. Then there is a constant Ci < oo 
with 

\\'^\\c°{Ti) — Ci||u||w/i.p(a)- (17) 

Furthermore we need the traces of functions from W^'^(O) on the boundary 90. 
Definition 4 
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Let O be a bounded domain with '^-boundary. That means the boundary 90 
is locally given by a function with a Lipschitz continuous derivative. Suppose p > 1 . 
Then the trace space of functions from W^'T'(fl) is defined as (see jH], pp- 37-38) 



W^^p '^(30) ={ue W^'P(9a) 



30 30 

for all a 



dSxdSy < oo 



[a-\, az) G Nq, oci + (X2 < 1}. (1^ 



If p > 1 , then the trace operator 

Tr : W^'^iCl] 
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(19) 



Tr : u 1-^ ulaa (20) 

is well defined in the classical sense. 

If the time t is fixed, then we consider the field variables to be the functions, 
which are defined in D. and belong to the proper function spaces. We will use the 
following abbreviations of function spaces and subsets: 

• Xp := (Lp(Q))^ for the volumetric loads, 

• Yp := W^'P(Q) for the components of creep strain tensor, 

• Yp := (Yp)"^ for the creep strain tensors, 

• Vp := (W^'P(Q))^ for the displacement fields, 

• Vp := (W^''P(n) n Wq'^(O))^ for the displacement fields with a vanishing 
boundary values (which correspond to the solid clamped at the boundary), 

• yPi'P^ ■={(jo eYp -.0 < (Jo[x) < 1 - |3i, IIcuIIyp < Pi} for the damage fields, 
where > (3i > 0, |32 > are fixed constants. Accordingly to Corollary 3.1, 
uj[x) is well defined and Y^' is a closed subset of Yp. 

Remark. The condition < uj(x) < 1 — |3i is natural to guarantee that the 
elasticity tensor © is positive definite. The second condition ||cu||yp < Pi imposes 
additional constraints both on the damage field and on the damage gradient. 



3.2 Reduction to zero prescribed displacements 

In this subsection we reduce the boundary value problem ([T)) — ffTUj) to the case of 
zero prescribed displacements along the boundary 9Q. 
Theorem 3.2 

Suppose that displacements, which are given on the boundary, satisfy the follow- 
ing smoothness condition: u* G (W^^p '^(90))^. Then there is a function H G Vp 
with aiaci = u* in the trace sense (see [8J. 

In what follows, we designate H by the same symbol as u*. 

Now we reformulate our problem in a standard way. We search a vector 
(u, £'^^, cu) G Vp X Yp X Yp, such that (u + u*, e'^^, cu) G Vp x Y^ x Yp is a solu- 
tion of dH) — dfoi). 
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3.3 Compact form of evolution equations 

In this subsection we rewrite the evolution equations © in a compact form 

£""(x,t) =7^(p(u,e"^a))(x,t)) (21) 

cb(x,t) =5(p(u, e^^cu)(x,t)). (22) 
To this end we introduce for every (u, e'^^, cu) G Vp x x Y^' 

P(U,£^\CU) := (£ii(u),£22(u),£i2(u),£^;,£5^,£^^,a)) G YJ, (23) 

£(u) :=l(Vu + VuT) G Y4. (24) 

For every p G M'' we define 

7^(p) := s(p) (cTvm(p))"-^(1 - Pz)"", (25) 

5(p]:=B ((T.M(p])"^(^-P7)-^ (26) 

CTvm(p) := P(cTll(p), 0-22(p), CTl2(p]), (27) 

s(p):=cr(p)-ltr((T(p)]I, (28) 

-fp' - ( Z Z ) 'p'> 

0n \ £ / 1 ^ \ / pi - P4 

Oil Mp):=(1-P7)y^^ ^1 P2-P5 |, (30) 

(Tl2 / V ^ / V 2(P3-P6) 



P(zi, Z2, Z3) := yzf + — ZiZi + 3z| for every zgM. (31) 

The constitutive relations ()30|) . (jHT|) are obtained from the general 3D relations 
under plane stress assumption (cti3 = O23, = Gss = 0). 

4 Main result 

The existence and uniqueness theorem states that a unique smooth solution to the 
initial boundary value problem (P) — (jlUj) exists in a certain time interval. 
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4.1 Formulation of the main theorem 



Theorem 4.1 

Let CI be a bounded domain with C^'^ -boundary, p > 2, T > 0, q G C°([0, T] , Xp), 
u* G Vp, G Yp, and cuq G Y^'-I^^. Then there exists J^ G (0,7] such that for any 
T' G (0,Ti] there is a uniquely determined mapping [u,t'^^,uj) G C'^([0, T'], V° x 



Yp X Yp) such that 



V- ((1 -cu)C(£(u + u*]-£n) =-q(t) for all tG[0,T'], (32) 



t 

W,Cv)[t) = (£S^CUo) + 





(7^(p(u + u*, t", cv),S[p[u + u*, £^^ o))) (s)ds (33) 



for every t G [0,T']. Here the evolution operators TZ,S are defined by (PT|) — ()3ip. 
Moreover, 

cu(x,t]<1 for all (x,t] G a X [0,T'], (34) 
(£-,a;]GCM[0,r],YjxYp], (35) 

£"^ = 7^(p(u+u^£^^cu), cb = 5(p(u+u^£^^cu), (36) 

£-(0) = £S^ a;(0) = a;o, (37) 

||£^-(t)-£S-||y4 + ||a)(t]-a)o||Yp <min( ^^^^^^^^ ,^) VtG [0,r]. (38) 

Theorem 4-1 is proved in the next section. 
Corollary 4.1 

// t/ie so/zc? geometry, applied loads, prescribed displacements, and initial data are 
smooth, then the fully coupled creep-damage model predicts a nonzero lifetime t* of 
the structure with a lower estimate Ti from Theorem 4-i (t* > T-\). 

Remark. Theorem 4-i assures that t* > Ti > only for smooth domains without 
notches. Numerous examples of FEM simulation of notched specimens show that 
the predicted crack initiation time t* tends to zero as the mesh-size decreases ([18 ). 

Remark. The lifetime estimate Ti depends on the constants |3i, Pi. Moreover, 
as it will be clear from the proof of Theorem Ti = Ti((3i/|32). It is natural 
that Ti — > as |3i — > since the lifetime of the structure made of almost broken 
material (mi_n(1 — Wq) — > 0) is negligibly small. Furthermore, Ti tends also to zero 

as pi tends to infinity even if |3i is finite. The physical interpretation of this result 
could be the following. The rupture time can be negligibly small in the case of 
big gradients of damage (||cuo||yp oo) even if the initial damage itself was not 
substantial (min(l — a>o) ~ 1). 

Example is provided in the subsections 4.2 showing that the dependence of t* 
on (32 can be interpreted as lifetime reduction due to damage localization. 
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4.2 Counterexample: lifetime reduction due to local imper- 
fections 

Consider a solid loaded by prescribed displacements on it's boundary as shown on 
figure HJ Assume that the boundary and prescribed displacements are smooth. We 
set the creep strain rate to zero (B = in 0). Consider a curve I of length L within 
the solid. Suppose that initial damage is concentrated near the curve I (see fig. Hj) 
and the initial creep is zero 

,^ h, — distfx, I] , ri- ^ , s 

a)o(x) := max(0, e^^ := 0. (39) 

It is obvious that 

min(l — uJo) = 1/2, HcdqIIyp ^ oo as H ^ 0. (40) 
We assert that 

t* ^ as H ^ 0. (41) 

To prove this assertion we can use the same argumentation as used in ^H]. The 
main reason the lifetime is decreasing is because the stress concentration factor near 
the curve tip tends to infinity as h. — > 0. 

5 Proof of Theorem ^.1 

To prove Theorem 4-1 v^e need several lemmas. 



5.1 Equilibrium equations with respect to u with a given 
and CD 



Define the family {Cw} P'2p, operators of linear elasticity 

C^:Y^^ Xp (42) 
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by the rule 

£^(u):=V-((l-a))Ce(u)) (43) 

for w e Y^'^^\ UG Vp. 
Lemma 5.1 

JLl 2P2 

The operator is bounded for all w G Yp^ ' . Moreover, the problem 

CM = -q (44) 

— 2(32 

has a unique solution u G Vp for all w G Yp^ ' ^ , q G Xp and 

||u||vp < C5,||q||xp. (45) 
Here C5 i < 00 does not depend on q. 

Proof. The boundness of C^^ follows from the following computations 

||/:a,(u)||xp=||V-((l-a;)C£(u))||x, 

< ||Vcu||xp • ||C£(u)||co + ||1 -culico ■ ||V- (Ce(u))||xp 

< C(|3i,(32)||u||vp. (46) 

Up to the rest of the article the expression Qi < C • Q2 should be understood as 
follows. The quantities Qi and Q2 are related to each other in such a way that there 
is a suitable constant C < 00, which depends only on [CI, E,"v,p,n, m, q, (3i, ^2) and 
Qi < C-Q2. 

For the proof of solvability of (jUj) and for estimate see ^ (p. 241). Partic- 
ulary, we have the following inequality 

||u||vp <C(a,E,T,p)^||q||xp. (47) 

Pi 

The lemma is proved ■ 

Remark. The influence of the damage localization on the strain localization is 
taken into account by (jTTj) . Indeed, 

C5 1 — > 00, as > 00. (48) 

Pi 

We denote by Cw\v° the restriction of Cw to Vp. Let be the inverse to Cuj\vo . 
Since fl^3j) holds, we see that 

ll^c.l<C5.i. (49) 

Lemma 5.2 

There is a constant C5.2 such that 

||cui - a)2||Yp, (50) 

ll>Ca,l -CZ:]\\ < C5.2IICU1 - CU2||Yp, (51) 
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\c° 



Pi 262 

for all cui, CU2 G Yp-^ ' 
Proof. Obviously, 

\\C^^ - Cao2\\ = sup ||V ■ ((a)2- a)i)Ce(u))||xp 

ll^llvp =1 

< sup (||V(a;i-a;2)||xp-||C£(u] 

l|u||vp=1 ^ 

+ llcul -CU2IIC0 • ||V- (C£(u))||xp) < C llcul -a)2||Yp. (52) 

Furthermore, 

ll^cul -^culll < lldl • ll^culll • ll^o. -C^A < C llcUl -CU2||Yp. (53) 

The lemma is proved ■ 
Lemma 5.3 

Let u* G Vp, T>0, q G C°([0,T],Xp], e^^ G C°([0, T], Y^), 

— 2(32 

cu G C°([0,T], Yp^ ' ). Then there exists a uniquely determined mapping 
U = U(u^£=^cu,q) G C°([0,T],V°) such that 

/:^(t)(U(t)) = -/:a,(t)(u*) - q(t) + V ■ ((1 - a)(t))Ce^^(t)) (54) 

for t G [0,T]. Moreover, 

||U(u^£^^a),q)||vp,oo < C5.3(||u*||vp + ||£'1y4,oo + Nllxp.oo), (55) 

l|U(u*, e-^ , cu\ q V U(u*, cu^ q^) Hvp.oo 

< C5.3(||CU^ -CU^||yp,oo(||u*||vp + W^\\y^^,oo + Nlxp.oo) 

+ ||£-^-£-^||Y4oo + ||q^-q^||xp,oo). (56) 



Proof. We claim that the mapping U G C°([0,T], V°) is uniquely defined by 
Indeed, at each instant of time we have by Lemma 5. 1 

U(u*, cu, q) (t) = ( - /:.,(t)(u*) - q(t) + V ■ ((1 - cu(t))Ce-(t))) . (57) 

Estimate ()55|) follows from 

Combining (jl^ . (jSnj, and (jSH) we note that 

\\C-^,C^^u*-C-J,C^^u*\\ < C llu^llvpllcu^-cu^Yp. (58) 

Note also that 

\\C-Jg' - C-Jg^ < C\\co' - cu^||Yp||g^||xp + C||g^ - g^lUp, Vg\ g Xp. (59) 



Substituting (^-qHt)+V - ((1 - cuHt))C£'="(t) j for gMn (jSHl) and combining §^ 
with (jSHj) . we get This completes the proof of Lemma 5.3 M 
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5.2 Evolution of cu and e'^^ with a given u 

Let us first analyze evolution operators IZ, S, which are defined by (j2SI), (j^ - 
Lemma 5.4 

There is a constant C5.4 with the properties to follow. For all i, j G {1 , 2}, 
k G {0, 1 , 7}, p\ G such that p}, P7 < 1 — ■'y we have the following estimates 

|D„7^tJ(p^)| + |D^5(p^)| < C5.4(1 + |pi|)--(-'-\ (60) 

|D^7^tJ(p^) -Dk7^,J(p^)| + \D^S[p')-D^S[p^)\ 

< CsaO + Ip^I + |p2|)'^^-("^.^)|pi - p2|. (61) 

Here D]^ = for k G {1 ,2, 7} and Do is the identity operator. 
Proof. First let us note that for all q > 2 

|P^(z)| < C|z|^ IP^tzj-P^lz']! < C(|zl + |z'l)^-V-z'|, (62) 



|V,P^(z)| < C|z|^-\ IV.P^tz)- V,P^(z')| < Cdzl + lzD^-^lz-z'!, (63) 

where P(z] is defined by (jSD). 

We also note that for all p\ p^ G R^, i, j G (1 , 2}, k G {1 , 2, 7} 

|cTij(p)| < CIpI, |0t,j(p^) - (Ty(p2)| < C(1 + Ip^l + Ip^Dlpi - p2|, (64) 

l^cTij(p)l < C(l + Ipl], l^(at,j(p^) - ay(p2))| < C|p^ - p^l. (65) 
9pk 9pk 

Here crtj(p) is defined by (jHUj) . 

The lemma is proved after some simple computations. Let us prove for example 
that 

|Di,5(p^) - D^S[p^]\ < C (1 + Ip^l + Ip^l)"^ Ip^ - p\ (66) 
We remark that S[p) has the form 

5(p)=P"^(ff,j(p])F(p7) (67) 

with F G C°°[0, 1 - Here ffy = [a^,^,az,2, (y^,2V■ 

|Dk5(pM-D^5(p2]| 

< |D4P-(a,j(p^))F(p}) - P"^(cTy(p^))F(p^))| 
+ |D^(P-(ay(p^))F(p}) - P-(cT,;(p^))F(p^))| 
< C |D^(P-(ay(p^)) - p-(ay(p^)))| + C |P"^(ay(p^)) - P"^(ay(p'))| 
+ C |D^P™(cTy(p2))| ■ |p^ - p^l + C |P"^(at^(p'))| ■ |p^ - p'l 

^ A + B+ C+ D. (68) 

We abbreviate 

Di,ct(p) :=DJcti,i(p),(T2,2(p),cti,2(p))"'. (69) 
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Further, 

A = C |V,P™(cTy(p^))D^CT(pi) - V,P"^(fftj(p'))Dfca(p2)| 

< C |V.(P-(ay(p^)) - P"^(cT,^(p')))| • |D^cT(pi)| 
+ C|V.P"^(cT,j(p^))|-|D^(a(pi)-a(p^))| 

< C (lay(p^)l + lay(p^)l)"'"' ■ lay(p^) - atj(p^)| ■ ID^alp^)! 

+ C|ay(p2)r-^|D4cT(p^)-cT(p^))| 

<c(i + ip^i + ip^np^-p^i. (70) 

In the same way we obtain 

max(5, C,D]<C[^+ |p^| + Ip^]]"" |p^ - p^]. (71) 

Combining this with (|7Up we get . The lemma is proved ■ 
Lemma 5.5 

There is a constant C5 5 with the following properties. For all M, T > 0, u\ G 
C°([0,T],Vp) with ||u^||vp,oo < Mfor lG{1,2}; e^^\ e^^^ G C°([0, T], Y^); oj\vj^ e 
C°([0,T],Yp) such that 

a^(x,t)<l-i^ V(x,t)Gax[0,T],lG{l,2} (72) 

||£-^(t) - £^^H0)||y4 + \\uo\t) - 0)^(0) II Y, < 1 V t G [0,T],l G {1,2}. (73) 
We abbreviate (recall (f^ ) 

p^:=p(u\£-\cu^ (74) 

Then 

7^(p^)(t)GYJ, 5(p^](t]GYp, VtG [0,T],IG{1,2}, (75) 

l|7^(p')||Y4,oo + l|5(p^)||Y„oo < C5.5(M+||£-^(0)||y4 + ||cuH0)||y,+1)"'^'"''"'^' (76) 
for I G {1 , 2}, and 

||7^(p^) -7^(p^)||Y4,oo + ||5(p^) -5(p^)||y„oo 

2 

< C5.5(m + }^ (||e-H0)||Y4 + W^'mW) + 1^ 
1=1 

|u^ -u^||vp,oo + ||e^^^ - e^^iY4 00 + W - c^1Iyp,oo ). (77) 



max(rn.,rL)+l 



Proof. This lemma is proved by Corollary 3.1 and Lemma 5.4- Let us estimate, 
for example, the value ||K||]_p(a)(t), where 



K(x,t) := 



g^(7^tJ(p^)-7^t,j(p^))(x,t]|, (7^ 
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with r, i, j e {1 , 2}, (x, t) e O x [0, T]. We obtain 



K(x, t)<Y_ |D.7^y {pUx, t]) _ D,7^,, (p^(x, t)) 

lc=1 ^ ^ 



< Y_ |D^7e,,^(pi(x,t]) -Dk7e,,j(p2(x,t)) 



ap[(x,t) 



9xt 



+ _^|D^,7^,,^(p2(x,t)) 



k=1 



9p[(x,t) dp^(x,t] 



dXr 9Xr 

7 

<^Ak-Bk+Ck-Di 



lc=1 

But by Lemma 5.4 and by Corollary 3.1 we get 

Aic < C (1 + |pi (x, t] I + |p^(x, t) I pi (x, t) - p^(x, t) I 

< C (1 + ||p^||y7,oo + ||p'(x,t)||YZ,^)'"^'^'"V^ - p'IIy^,oo 
We get by the same argument 

C^< C |D^7^y(p^(x,t))| < C (1 + ||p^||y^,oo)"^^'^'"^ 

Evidently, 



llBiclk, = 



ap[(x,t) 



iDkllu = 



9Xr 

ap^(x,t) ap^(x,t] 



axr 



ax. 



Lp < IIp^ -p^||y7,o 



Hence 



k=1 



2|i \ max(Ta,n)+1 m i 2\\ 

• IIP - P ||y7,oo 



k=1 

<C (1+||p^||y7,oo + ||p'||>v,oo) 
It remains to check that 

IIP^ - P'||>-,oo < C (||U^ -U^||v,,oo + lle'^^^ - £'''||y4,oo + W - a^'||Yp,oo) 

and 

+ ||p^||y^.oo + ||p'||y^,oo) < C (m + ^ (lle-HO)llY^ + ||a)HO)||Yp) + 1). 

1=1 
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This concludes the proof of Lemma 5.5M 

We now prove the existence of the solution (£'^^,cu) G C°([0,T],Yp x Yp) if the 
displacement field u G C°([0, T],Vp) is given. 

Lemma 5.6 

Let M > 0, G Y^, and Wq G Y^i -P^. Put 



C5.5(M + 2||.§1,. +2||cuo||v, +l)"'^-'"^"'"^2 . (87) 
^ ^ min((3i,2|32j-l 



with C5.5 from Lemma 5.5 and Ci from Corollary 3.1. Let T' G (0,To] and 
u G C°([0, T'], Vp), such that ||u||vp,oo < M. Then there exists a unique mapping 
(e^^a)) = (£^^cu)(u, £S^a)o) G C°('[0, T'], Y^ x Yp) such that 



(7^(p(u,£^^a))),5(p(u,£^^a))))(s)ds, (88) 



||£-(t) - £-(0)||y4 + ||a;(t) - a)(0)l|Y, < V t G [0,1']. (89) 

Moreover 

a;(t) G Yp^'""^ VtG[0,T'], (90) 



[£- a;]GCM[0,r],Y>Yp], (91) 

(e- d))(t) = (7^(p(u,£-,a))),5(p(u,£- a))))(t) VtG (0,7'). (92) 

Proof. As it is done in j2| for partly coupled damage model, we adapt the 
standard proof of the existence of solutions to ordinary differential equations in 
Banach spaces. We define the closed subset of C°([0,T'],Yp x Yp) by 



M := {[t'\cu) G C°[[OJ']X X Yp) : (£^a))(0) = [tl\cuo], 

li 

~2(1+Ci 



£-(t) - £-(0)||y4 + ||a)(t) - cvmW, < ""^l+'c^^"^ }- (93) 



The application of Corollary 3.1 yields for every (£'^^, cu) G TVf 

cu(t) G Y^'^*^' VtG[0,T']. (94) 
Note that for all t G [0, T'], (£=^ uj] e M 

||£-(t) - £-(0)||y4 + ||a)(t) - a)(0)||Y, < "^^if'^l^f < 1- (95) 
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Hence, by Lemma 5.5, we obtain for t < t', (e'^^, cu) G 
t' 

|(7^(p(u, cu)),5(p(u, £-,a;)))(s)||Y4,y/s 



< C5.5(M+ ||£-^(0)||v. + ||cuHO)||y, + l)'^^^'"^'"'^\t-t' 



min(|3 1,2(32) , . 



Let the mapping T : ^ C°([0, T'], Yp x Yp) be given by 
T(£^^a))(t) 



t 



;£S^^o) + 



(7^(p(u,£^a;)),5(p(u,£^^cu)))(s)ds VtG[0,r]. (97) 



Accordingly to (j96j) . this mapping is well defined. Taking into account (j96|) . we 
obtain 

||r(£^^a))(t)-(£S^a)o)||y4>,y^ < "^'^^[^I'^^f VtG[0,r]. (98) 

Therefore, T[Ai) C TW. By Lemma 5.5, it follows that T is a contraction with 
respect to the norm of the space C°([0,T'],Yp x Yp). Indeed, for every instant of 
time t we have 



||r(£-\a;^)(t)-T(£-^a;2)(t)||Y4,Y, 

2 



/ , — , , , X \ max(m,TL)+1 

< c5.5(m + Y_ (i|£^'Ho)iiy4 + wmw^) + 1) 
1=1 

• (||£-^ - £^^'||y4,oo + IIo^' - a^'l|Yp,oo)r < 

1(||£-^ - £-iY4,oo + IIO^^ - a;'||Yp,oo). (99) 

It follows from Banach's fixed point theorem (see, for example, jS]) that there exists 
a unique (£'^^, w] E Ai such that T{t'^'^, uj] = [t'^^, uj). Thus, we have proved that 
the pair (£^^ w] G C°([0,T'], Y^ x Yp) is uniquely defined by pl). 

To conclude the proof it remains to note that the mapping t h-> 
(7^(p(u, £'=^, cu)),iS(p(u, £'^^, cu)))(t) is continuous from [0,T'] to Yp x Yp. Since 
fl90|) holds, we may use Lemma 5.4 to prove that 

l|(7^(p(u, £- a))),5(p(u, E-\cu)]){V] 

-(7^(p(u,£-,a;)),5(p(u,£- a))))(t2)||Y4xYp ^0, as ^ t^. (100) 

Lemma 5.6 is proved ■ 
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Now we need to estimate the difference between two solutions of (|8a|). (|89|). Let 
us abbreviate 

e = C5.5(M + 2||e-||Y4 +2||a;o||Y, + 1)°-^'---'+^ (101) 

with C5.5 from Lemma 5.5 and Cj from Corollary 3.1. 
Lemma 5.7 

Let M > 0, K > 1 , G Y4, and cuq e Y^' -P^. Put 



T 



1 •= 



e-2max(K,— (102) 
mm((3i,2|32j J 



where t is given by (ITTITll . Let T' G (0,Ti], and u\\i^ e C°([0, T'], Vp), such that 
||u.^||vp,oo < M for I G {1,2}. Assume that 

(£-\a)^) = (e- a))(u\£- cuo] G C°([0,T'],Y4 x Yp),l g {1,2} (103) 

are defined by dHHI), dHHI). Then 

II e^'-' - £^^'||y4.oo + IIo^' - cu^llYp.00 < ^||ui -u^llv^oo. (104) 
Proof. By Lemma 5.5 we have 

||7^(p^) -7^(p^)||Y4,oo + l|5(p^) -5(p2)||Y„oo 

< e • (||u^ -u^llvp.oo + lie"'' - e''iY4,oo + l|a>^ - cu^Hyp.oo). (105) 
Therefore, since ()88|1 holds, we obtain 



I fCrl £.ci'2|| I / 1,2 1 1 

|£ — £ IIy4,oo + 11"^ — llYp.oo 



< eT'(||ui -u^||vp,oo + ||£'^^ - £''1y4,oo + llo^^ - a^'||Yp,o 

- 2!<"''' """'"^-^ + ^(ll^"' - £"'IIy4,oo + W- C^'IIYp.oo). (106) 
Now inequality ()104j) follows from ()106|) . Lemma 5. 7 is proved ■ 

5.3 Equilibrium equations coupled with evolution law 

In this subsection we solve problem (jS2I), ()33|1 and prove Theorem 

Proof of Theorem 4.1. We choose Ti as in Lemma 5.1 where we put 

K := 2 ■ C5.3(||u*||vp + ||£Slv4 + ||q||xp,oo + l), (107) 

M := C5.3(||u*||vp + ||£S'llv4 + ||q||xp,oo + l) + ||u*||vp (108) 

with C5.3 from Lemma 5. 3. In order to use Banach 's fixed point theorem we define 
the closed subset of C°([0, T'], V°) by 

M:=[ue C°([0,T'],V°) : ||u + u*||vp.oo < m}. (109) 
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Let the mapping T : M ^ C°([0, T'], V") be given by 

r(u) = U(u*, £^^(u + u*, cuo), a;(u + u*, cuo), q) , (HO) 

where (e'^^, cu) (u, t'^ , cuq) is defined by (jHHI) and U(u*, e'^^, cu, q) is introduced in 
Lemma 5.3. 

Let us show that T[}A] C C°([0, T'], Vp]. In fact, accordingly to Lemma 5.6, 
for all u e M we have (e^^ cu) (u + u*, t^^cvo) E C°[[OJ'],Y^ x Yp). Therefore, 
the assertion follows from Lemma 5.3 and the mapping T is well defined. 

Now let us show that T[M] C M. Using (j^SI) and (jEHl) we obtain for u G 

||T(U) +U*||vp,oo < ||T(u)||vp,oo + ||u*||vp 

< C5.3(||u*||vp + ||e^^-(U + U*, eS'.^0)||Y4,oo + Nllxp.oo) + ||u*||vp < M. (Ill) 

Let us prove that T is a contraction. Taking into account ()1()4|1 . and the choice 
of K, we obtain 

||T(ui)-T(u2)||vp,oo < 

C5.3(||CU^ - Cu2||Yp,oo(||u*||vp + ||£o'IIy;; + l|q||Xp,oo + l) 
1 II _ct1 c.ct2|| \ 

< %(N*llvp + ||£o'IIy4 + ||q||xp,oo + l)||u^ -U^llvp.oo 

< l||u^ -u^llvp.oo. (112) 

From Banach's fixed point theorem it follows that there is a uniquely determined 
mapping u G A4, such that T[u) = u. To obtain the mapping (u, £'^^,cu) G 
C°[[0,J'],V° X Y4 X Yp) that satisfy §^ and (jMl), we put (see Lemma 5.6) 

£^^:= £^^(u + u^£^^a)o), a):=a)(u + u^£^^a)o). (113) 

The existence of the solution to (jH^ . (jH^ is proved. Although the uniqueness in Jvi 
is guarantied by Banach 's fixed point theorem, it remains to check that the solution 
is uniquely determined by Assume the converse, then there are two 

different solutions of (jSH), (jH^ 

(u\£-\OGC°([0,T'],Vp°xY4xYp), IG{1,2}. (114) 

Put 

to:=max{t G [0,7'] : (u\ £^^\ cu^)(t) = (u^ £^"^ cu^) (t) for t G [0,t]}. (115) 
Hence, 

(u\£^^\cu^)(t) = (u^£^^^cu^)(t) fortG[0,to]. (116) 
Moreover, to < T' and for every t G (to,T'] there exists t G (to,t], such that 

(u^£-^cul)(t)^(u^£-^a;2)(t). (n?) 
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Arguing as above, we prove the uniqueness of the solution 
(u, e^^ tu) e C°([to,T"^^'], V° x x Yp) to the new problem 

V- ((1 -cu)C(£(u + u*)-en) =-q(t) V t G [to, T""'], (118) 
(£- a))(t) = (e-^(to),cunto)) 

(7^(p(u + u^e'^^cu),5(p(u + u^e^^cu))(s)ds Vte[to,r^^'] (119) 



t 

+ 

to 



||u(t] +u*||vp < M"'^" Vt G [tcT"'^"'], (120) 
with some new parameters jnew/^j>^new rjj^g reader will easily prove that there 
exists t e (to.minlT^*^"',!')], such that 

l|uHt)+u*i|v^ < M"*^™ VtG [to,t], lG{1,2}. (121) 

That means, that 

(u\e^^\a)^)(t] = (u^£^^^a)^)(t) Vte[to,t]. (122) 
This contradiction proves the theorem ■ 

6 Conclusions 

The creep damage problem is formulated in a well-poscd manner. Theorem 4-i 
states that a imique smooth sohition to the Kachanov-Rabotnov problem exists in a 
certain time interval [0, TiJ. The corresponding function spaces Xp, Yp, and Vp reflect 
the essence of the system of equations and can be used for a proper mathematical 
analysis of the problem. Particulary, clear definitions of terms "stable", "unstable", 
and " convergency" can be given. 

It is shown that the requirements of the existence theorem {Theorem 4-i) have a 
physical meaning and the violation of these requirements directly affects the lifetime 
estimate. 

If we do not impose any restrictions on the gradient of initial damage (such 
situation corresponds to (Si = oo), then the lifetime t* of the structure can be 
arbitrary small even if min(1 — tUo) > |3i > 0. 

This damage localization effect is characterized at each instant of time by the 
quantity 

A(t) = J^^^^lt). (123) 
mm( I — w) 

The value A(t) controls the remaining hfe of the structure trest '•= t* — t. 

trest — > 0, as A— >oo. (124) 

Thus, the estimation of A gives an answer to the question when the damage becomes 
critical. This measure of damage localization can be adopted to improve monitor- 
ing and inspection strategies used to secure the reliable operation of engineering 
structures. 
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